The Moving Glass Phase of Dirty Type II Superconductors 
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We study numerically the motion of vortices in dirty type II superconductors. In two dimensions 
at strong driving currents, vortices form highly correlated "static channels". The static structure 
factor exhibits convincing scaling behaviour, demonstrating quasi long range translational order in 
the transverse direction. However order in the longitudinal direction is only short range. We clearly 
establish the existence of a finite transverse critical current, suggesting strong barriers against small 
transverse driving forces. These results are most consistently interpreted in terms of the recently 
proposed moving glass picture, modified to account for the strong anisotropy. 
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Many condensed matter systems reach higher levels of 
organization by forming periodic media. Examples range 
from crystalline solids to Wigner crystals, charge density 
waves and vortex lattices in type II superconductors. A 
central issue is the effect of disorder on the stability of 
such systems. In the case of flux lattices, for weak disor- 
der it seems sufficient to concentrate on elastic deforma- 
tions which were argued to lead to a power law decay 
of lattice correlations, when the periodicity is properly 
taken into account. This phase was christened a "Bragg- 
glass" [jjj. The irrelevancy of topological excitations has 
been recently confirmed in three dimension, while the 
situation is marginal in two dimension g]. For strong 
disorder topological excitations such as vortex loops can 
become relevant, giving rise to a Vortex glass ||. While 
numerous additional scenarios have been proposed [Q, 
there is growing experimental evidence supporting the 
basic picture of two types of glasses as the disorder or 
the magnetic field is increased |]||. 

Upon increasing the external force beyond a certain 
depinning strength, these periodic (or glassy) media be- 
come mobile. Early work developed perturbation stud- 
ies at high velocities v in powers of l/v . Recently 
Koshelev and Vinokur argued that the effect of the ran- 
dom potential is seriously weakened at high velocities, as 
it "averages out" , thus reestablishing the long range solid 
order || . Their numerical simulations suggested that this 
change occurs abruptly, giving rise to a genuine dynamic 
phase transition. On the other hand Giamarchi and Le 
Doussal pointed out that some components of the disor- 
der remain unsuppressed, as they represent static pertur- 
bations pOj . These destroy the moving solid, and stabi- 
lize a glassy phase instead with quasi-long range order 
(QLRO) only, giving rise to a moving glass. Balents and 
Fisher, in their study of the related charge density wave 
systems also find that the moving phase possesses QLRO 
only, but may nevertheless support true long range tem- 
poral order |ll[] . The physics of the moving glass phase is 
|Tof that the vortices form highly correlated static chan- 
nels. This picture leads to a logarithmic decay at large 



distances for the displacement correlations. It also gives 
rise to diverging potential barriers and consequently a fi- 
nite critical current against an additional transverse cur- 
rent. 

In this Letter, we report a detailed numerical study 
of the moving glass phase. We establish that the ba- 
sic picture of "static channels" is indeed characteristic 
of the system. By studying the static structure factor 
we demonstrate the existence of quasi long range trans- 
lational order with an algebraic decay in the transverse 
direction but find only short range longitudinal order, 
giving rise to a very anisotropic glass. As a direct con- 
sequence we predict the absence of a narrow band noise 
for moving glasses in the 2D vortex systems. Finally the 
existence of a critical transverse current will be clearly 
established. 

We employ overdamped Molecular Dynamics (MD) 
simulations at zero temperature to study two dimensional 
interacting vortices in the presence of point disorder, 



(i) 



Here 7 is the damping parameter, Rj specifies the pin- 
ning center positions, i\ denotes the location of the i-th 
vortex, and Fl is the Lorentz force, exerted by the exter- 
nal driving current. The force between vortices is given 
by 



F„(r) - F (l - r 2 ) 2 



(2) 



where f = r/R cut ,F = V /R cut , and we choose R cut = 
3.6ao, where ao is the mean vortex spacing. Here 7, do, 
and Vo(= < &q/87t 2 A 2 ) define the units of time, length and 
energy, respectively. The pinning force is taken as 



t (f) = _4F p (l-f 2 ) r 



(3) 



Here r = r/R V i n and R P i n — 0.25ao. We worked with 
a fixed density of pinning centers of 5.77/a 2 ,. This set 
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of parameters was chosen to optimize the convergence of 
the numerical procedure. 

Now we construct the phase diagram in the driving 
force - pinning strength plane, at zero temperature. With 
increasing driving currents three phases emerge: a pinned 
glass, a plastic flow regime, and some kind of an ordered 
phase. At low forces the vortices remain pinned, form- 
ing a glassy phase. As the Lorentz force is increased 
beyond a critical value Fd, the vortices depin. This tran- 
sition, and in particular the value of Fd can be well cap- 
tured by studying the current-voltage (IV) characteris- 
tics. The resulting values of Fd were used to construct 
the lower phase boundary in Fig.l. For strong disorder, 
just above Fd vortices form a pattern of pinned and un- 
pinned regions, often described as "plastic flow" 0. In 
this regime Fd scales linearly with the pinning strength, 
whereas for weak disorder the relation is quadratic. The 
near-linearity of the phase boundary in Fig.l indicates 
that we concentrated on the regime of strong disorder. 
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FIG. 1. The dynamic phase diagram. F p is the pinning 
strength, Fl is the Lorentz force. Open circles represent F g , 
solid circles Fd. 

Upon further increase of the driving force, Koshelev 
and Vinokur argued |J that the effects of the disorder 
on the lattice displacements "average out", leading to 
solid ordering in 3D, and quasi long range translational 
order in 2D. Alternatively, Giamarchi and Le Doussal 
suggested fl(| that the vortex system forms a moving 
glass. To distinguish between these propositions we first 
determine the location of the phase boundary, then ex- 
plore the physics of the high- velocity phase. 

The phase boundary F g between the plastic flow 
regime and the putative moving solid can be established 
by measuring the static structure factor SQs). In the 
plastic flow regime the absence of ordering manifests it- 



self in a central peak and a structureless ring (see lower 
panel in Fig. 2). In the high velocity phase Fl > F g , 
one expects to see six-fold coordinated Bragg peaks, if a 
moving solid is formed. We do indeed observe a sharp 
transition into a phase with well developed peaks, how- 
ever the peak pattern is strongly anisotropic (upper panel 
in Fig.2). 




FIG. 2. The static structure factor 5(k). F L = 0.6 in the 
upper and 0.2 in the lower panel. The number of vortices is 
960 and the disorder strength F p = 0.16. The critical force 
F g indicating the transition from the plastic flow regime to 
the moving glass is about 0.35. 

The upper phase boundary in Fig.l was determined 
by mapping out F g for several disorder strengths. One 
can see that for strong disorder indeed all three expected 
phases are observed, whereas for weak disorder there is 
no robust evidence for an intervening plastic flow regime. 
Either that phase occupies a very slim region in param- 
eter space, or there is a direct pinned Bragg glass-to- 
moving glass transition. This transition is much harder 
to identify because both phases exhibit quasi long range 
order, and thus the structure factors are very similar in 
the two phases. 
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The central issue of our paper is to elucidate whether 
at high velocities the system supports a moving solid or a 
moving glass. To address this issue we first analyse 5(k), 



S(k) 
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The pinning strength F p is fixed to be 0.16 and the ap- 
plied force Fl — 0.6 is well above the corresponding 
critical force F g = 0.35. We simulate five different sys- 
tem sizes with fixed vortex density and number of vor- 
tices ranging from 240 to 1500. The initial configurations 
are chosen randomly. We let the MD simulations evolve 
with time, make sure that the system reaches its steady 
state, then freeze the vortex configuration and measure 
S(k). In the steady state the vortices form an orderly 
array. Its principal lattice vector in most cases is aligned 
with the direction of motion. It was argued that the sys- 
tem chooses such an orientation to minimize the power- 
dissipation ||]. However the details of this alignment 
have yet to be understood. The peaks at the reciprocal 
lattice vectors perpendicular to the motion exhibit quite 
convincing power law behavior. The height of the other 
peaks however decay very rapidly with increasing system 
sizes, suggesting that the corresponding correlations are 
short ranged. 



with Sk = |k— Go|. In Fig.(||) the scaling function G(x) is 
plotted with respect to the dimensionless scaling variable 
x = 5kL for the five system sizes. The peak amplitude 
scales with the system size as l?~ Va with v s = 0.53 ±0.1, 
as shown in the inset of Fig. (||) . Using this value of v s the 
normalized structure factor exhibits a convincing data 
collapse onto a single curve. This confirms the scaling 
behaviour 5(k) = |k — G | -1 ' 47 around the peaks. 

This implies that phase slips in the transverse direction 
are strongly suppressed, justifying the elastic approach. 
In contrast, the peaks at momenta with nonzero longitu- 
dinal components decay rapidly with system size. Several 
physical mechanisms can lead to such a decay. Exten- 
sive study of single snapshots of the spatial distribution 
of vortices suggests that phase slips between longitudi- 
nal boundaries of elastic domains are primarily respon- 
sible. This observation questions the elastic theory for 
this direction. Clearly a more complete understanding is 
needed on this issue, especially in two dimension |T^| . 

A measurable consequence of the absence of transla- 
tional order in the longitudinal direction should be the 
corresponding absence of narrow band noise in 2D driven 
vortex systems. The same conclusion was reached for the 



analogous CDW models in 2D in Ref. |11 




FIG. 3. Finite size scaling of the static structure factor 
S(k). The driving force Fl — 0.6 is well above the critical 
force F g ~ 0.35. The inset shows power law dependence of 
the peak heights with varying system size; S P (L) ~ L 2 ~ Vs 
with u s ~ 0.53 ± 0.1. 

To study the peaks of 5(k) at the reciprocal lattice 
vector Go = (0, ±47r/v3)) we write down the following 
finite size scaling form: 

S(6k,L) = L d -^G(SkL) 7 (5) 




FIG. 4. The static channels in the steady state. The small 
fluctuations transverse to the channels are due to the pertur- 
bation from the rapidly varying time-dependent component 
of the disorder, viewed in the rest frame of vortices. 

The scaling behaviour of S(k) does not distinguish be- 
tween a quasi-solid and a moving glass. The proposi- 
tion of the glassy phase rests on the argument that cer- 
tain components of the disorder do not average out, but 
present a static perturbation. If so, the moving vortices 
should form static channels, which do not change their 
shape with time. To study this we now map out the 
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trajectories of the vortices. Making sure that the flow 
reached its steady state, we take a large number of con- 
secutive snap shots, which are then displayed on top of 
each other. The resulting Fig.(^) clearly demonstrates 
the formation of static channels. The very existence of 
these time-independent channels lends strong credence to 
the moving glass picture. 

The moving glass picture also implies the existence of 
diverging barriers against small transverse currents, lead- 
ing to a finite transverse critical current. This critical 
current is the largest when one of the principal lattice 
vectors is parallel to the motion |l(| . 

We select 50 disorder realizations which lead to a 
steady state with one of its primitive lattice vectors par- 
allel to the direction of the velocity. After the system 
reaches the steady state, a small additional transverse 
force is applied and the transverse velocity v y is mea- 
sured. Fig.(||) clearly exhibits a finite critical transverse 
force F% ^ 0.006 ± 0.001. F$ is much smaller than the 
longitudinal critical force Fd — 0.1, since the effective 
disorder strength is inversely proportional to the velocity 
of the moving glass JTo| . Simulations for several different 
system sizes result in identical critical currents, indicat- 
ing that this threshold behaviour is not a finite size effect 
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FIG. 5. IV characteristics, showing transverse velocity v y 
as a function of transverse force F v 



ring shape in the plastic flow regime into a phase with an 
anisotropic peak pattern. The Bragg peaks suggest the 
existence of power law order in the transverse directions, 
but only short range longitudinal order. Thus we expect 
the absence of a narrow band noise. We demonstrated 
the formation of static channels, and found a finite criti- 
cal transverse current. These results are consistent with 
the suggestion that the driven vortex system in two di- 
mension forms an anisotropic moving glass. 
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Here F L = 0.6. The 
0.006 ± 0.001 is clear. 
The straight line represents a free flux flow response. 



existence of a finite critical force F% 



In sum, we explored the moving phase of vortex sys- 
tems at high velocities. We measured the IV character- 
istics to determine the critical depinning force Fd, estab- 
lishing the phase boundary between the pinned glass and 
the plastic flow regime. Next we studied the structure 
factor S*(k), which exhibited a sharp transition from its 
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